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Abstract Statistical properties of order-driven double-auction markets with Bid-Ask 
spread are investigated through the dynamical quantities such as response function. 
We first attempt to utilize the so-called Madhav an- Richards on- Roomans model (MRR 
for short) to simulate the stochastic process of the price-change in empirical data sets 
(say, EUR/JPY or USD/JPY exchange rates) in which the Bid- Ask spread fluctuates 
in time. We find that the MRR theory apparently fails to simulate so much as the 
qualitative behaviour ('non-monotonic' behaviour) of the response function R(l) (I de- 
notes the difference of times at which the response function is evaluated) calculated 
from the data. Especially, we confirm that the stochastic nature of the Bid- Ask spread 
causes apparent deviations from a linear relationship between the R(l) and the auto- 
correlation function C'(l), namely, R(l) oc —C(l). To make the microscopic model of 
double-auction markets having stochastic Bid- Ask spread, we use the minority game 
with a finite market history length and find numerically that appropriate extension of 
', the game shows quite similar behaviour of the response function to the empirical evi- 

dence. We also reveal that the minority game modeling with the adaptive ('annealed') 
look-up table reproduces the non-linear relationship R(l) oc —f(C(l)) {fix) stands for 
a non-linear function leading to 'A-shapes') more effectively than the fixed ('quenched') 
look-up table does. 
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1 Introduction 

How a specific trading mechanism effects on the price formation is one of the essential 
queries to understand the process and outcomes of exchanging assets under a given 
concrete rule. To investigate the issue, a lot of studies concerning the micro-structure 
of markets have been done in various research fields [T|. 

Recently, lots of on-line trading services on the internet were constructed by several 
major banks such as the Sony Bank [2]. As the result, one can gather a lot of trading 
data sets to investigate the statistical properties extensively. As such studies, several 
authors focused on the fact that the Sony Bank uses a trading system in which foreign 
currency exchange rates change according to a first-passage process (FPP) 
Automatic FOREX trading systems such as the Sony Bank are now popular in Japan 
where many investors use a scheme called carry trade by borrowing money in a currency 
with low interest rate and lending it in a currency offering higher interest rates. With 
these demands in mind, several studies have been done to investigate the stochastic 
process and made a model of it to reproduce the FPP in order to provide useful 
information for customers [6ll7ll8ll9lll0lllllll2U13] . 

The data sets of the Sony Bank rate [2] are composed of time index and trading 
rate at that time. As we explained, a huge number of market data are reduced to a 
small amount of it, namely, the number of the Sony bank rate is reduced by the first- 
passage process and unfortunately, the market rates behind the Sony Bank rate are 
not available for us. 

As well-known, there are several data sets whose price are determined by the so- 
called double- auction system. In the double-auction market, each trader (investor) posts 
his (or her) selling price or buying price in market order of a specific commodity with 
its volume to the market. Then, the market maker determines the minimum price of 
buying orders, what we call Ask, and the maximum price of selling orders, the so-called 
Bid at each trading time and discloses these prices to the public. Then, the difference 
between the Bid and the Ask is referred to as spread or Bid- Ask spread. In market rates 
available for traders (on the web for instance), there are two types of Bid- Ask spread, 
that is, 'constant' or 'distributed', and which type of spread is disclosed depends on 
the market makers (securities companies). 

Results of market making, especially, statistical properties of the Bid- Ask spread 
might have an impact on the market and several studies have been done to reveal the 
relationship between the properties of Bid- Ask spread and behaviour of the market 
1 , 14, 15 , 16 , 17, 18 . For instance, Madhavan, Richardson and Roomans 15 proposed 
a phenomenological model to explain the price dynamics of double-auction market 
in market order, however, their model is apparently limited to the case in which the 
Bid-Ask spread remains constant during the price dynamics. Therefore, much more 
extensive studies including empirical data analysis seem to be needed to investigate to 
what extent the model proposed by Madhavan, Richardson and Roomans can explain 
the behaviour of market with stochastic Bid- Ask spread through some relevant quantity. 

In this paper, we investigate statistical properties of double-auction markets with 
Bid-Ask spread through the dynamical quantities such as response function. We first 
attempt to examine the so-called Madhavan- Richardson- Roomans model (MRR for 
short) to simulate the stochastic process of the price-change in empirical data sets 
(say, EUR/JPY or USD/JPY exchange rates) in which the Bid- Ask spread fluctuates 
in time. We find that the MRR theory apparently does not simulate so much as the 
qualitative behaviour ('non-monotonic' behaviour) of the response function calculated 
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from the data sets. It is possible for us to show that a linear relationship R(l) oc 
—C(l) between auto-correlation function C(l) and response function R(l) holds for 
the MRR model. Namely, these two macroscopic quantities are related each other and 
the relationship should be explained from the microscopic view point as statistical 
physics provides the microscopic foundation of thermodynamics. Moreover, we find 
that the linear relationship R(l) oc — C(i) is apparently broken down in order-driven 
double-auction markets with fluctuating Bid- Ask spread. This fact tells us that on the 
analogy of physics, the phenomenological MRR model for the constant spread might 
be regarded as 'thermodynamics' which usually deals with the macroscopic quantities 
such as price, auto-correlation, response functions and the relationship between them. 
It does not need to consider the detail behaviour of microscopic ingredients such as 
traders. In this paper, we show that the phenomenological model is apparently limited 
and fails to reproduce the dynamical quantities C(l),R(l) efficiently. 

Hence, here we attempt to construct a kind of 'statistical mechanics' in finance, 
which provides a microscopic foundation of phenomenological theory such as the MRR 
model. For this end, we utilize the minority game with a finite market history length 
having the distributed Bid- Ask spread to reproduce similar behaviour of macroscopic 
dynamical quantities as the empirical evidence shows. 

In our minority game modeling, we first fix each decision component (buying: +1, 
selling: —1) in their look-up tables before playing the game (in this sense, the decision 
components are 'quenched variables' in the literature of disordered spin systems such as 
spin glasses [19]). We also consider the case in which a certain amount of traders update 
their decision components according to the macroscopic market history (they 'learn' 
from the behaviour of markets) so as to be categorized into two groups with a finite 
probability (in this sense, the components are now regarded as 'annealed variables'). 
Namely, at each round of the game, if the number of sellers is smaller/greater than 
that of buyers, a fraction of traders, what we call optimistic group/ pessimistic group, 
is more likely to rewrite their own decision components from — 1/+1 to +1/— 1. We 
find that the minority game modeling with the adaptive look-up table reproduces the 
non-linear relationship R(l) oc —f{C(l)) (f{x) stands for a non-linear function leading 
to 'A -shapes') more effectively than fixed (frozen) look-up table does. 

This paper is organized as follows. In the next section [2] we explain our data sets 
and their format. We investigate their statistical properties. In the next section [3] 
we evaluate two relevant quantities, namely, the auto-correlation function C(l) and 
the response function R(l), which are our key quantities to discuss the double- auction 
markets, for our data sets. We find that a linear relationship R(l) oc —C(l) holds for 
the data sets having a constant Bid- Ask spread, however, the relation is broken for the 
data with a stochastic spread. In section[3] we introduce the so-called MRR model as a 
phenomenological model and derive the C(l) and R(l). The difference between the MRR 
theory and the empirical evidence, the origin of the difference is discussed. In section 
[5] wc introduce and modify the minority game with a finite market history length 
and apply to explain the typical behaviour of the response function for the data with 
stochastic Bid- Ask spread. In section [SJ we reveal that the minority game modeling 
with the adaptive ('annealed') look-up table reproduces the non-linear relationship 
R(l) oc —f(C(l)) {f(x) stands for a non-linear function leading to 'A-shapes') more 
effectively than fixed ('quenched') look-up table does. In section]?] we comment on the 
possible extensions of our approach. The last section is summary. 
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2 Statistical properties of data sets 

In order to check the validity of our modeling of markets, we gathered data sets of 
double-auction markets from the web site http://www.metaquotes.net/ 20 by using 
the free software MetaTrader^. We shall explain the data format of the MataTrader 4- 
We used the script which is available on the web [20]. By using the script, the data 
sets are stored as the following format: 



2009/12/24 , 17 : 17 : 40 , 131 . 053 , 131 . 092 
2009/12/24 , 17 : 17 : 41 , 131 . 053 , 131 . 088 
2009/12/24 , 17: 17:41 , 131 . 052, 131 . 088 
2009/12/24 , 17: 17:43 , 131 . 048, 131 . 071 
2009/12/24 , 17: 17:44, 131 . 043 , 131 . 076 



From the far left column to the far right, transaction time (Year/Month/Day, hour:min:sec), 
Bid, Ask are shown. For instance, the first line denotes the Bid is 131.053 and Ask is 
131.092 on 24th December 2009 at 17:17:40. In this paper, we treat the data set written 
by the above format. Among data sets concerning various different financial assets, we 
shall use here specific four data sets, namely, USD/JPY exchange rates (23rd-28th Oc- 
tober 2009), EUR/JPY exchange rates (22nd-28th November 2009), NasdaqlOO (22nd- 
31st October 2009) and price of gold (28th-30th October 2009). Each data set contains 
10 5 -data points. 

In the conventional (standard) data for continuous-time double-auction markets, 
we usually use the data having transactions (buying or selling price and the transaction 
time) including the quote (Bid and Ask prices posted to the market with the time). 
However, unfortunately, the data set provided by the MetaTrader^ does not contain 
any information about the transaction. Namely, the 'Bid and Ask values' we mentioned 
above are the best selling price and the best buying price, and the time at which the 
transaction takes place. Therefore, the price itself is not available for these data sets. 

Hence, here we assume that the mid point mt of the Bid bt and Ask at at time t, that 
is, mt = (at + 6t)/2 is a sort of buying or selling price when the transaction is approved. 
Then, we shall define the sign of the 'return' of the mid points (the difference between 
successive mid points) as a Selling- Buying signal ej, namely, et = sgn(mt+i — mt). Of 
course, these definitions of 'prices' and the 'Selling-Buying signals' are different from 
the conventional one, however, we shall try to investigate the behaviour of the system 
having such a slightly different definition of the prices and signals in limited data sets. 

In this section, we first calculate the histogram of the return of the mid point Amt = 
mt+i — mt- The results are shown in Fig. [1] From these panels, we clearly find that the 
return of the mid point is distributed with 'heavy tails' as the conventional return of the 
price has |21] , To compare the results with Gaussians, we calculate the empirical mean 
Am = (1/T) Ylt=o Amt and the empirical variance <7^ m = (1/r) Y]J—q (Amt — Am) 2 
for each data and plot the corresponding Gaussian J\f(Am, ~0Am) m the same panel. 

We next focus on the Bid- Ask spread St = at — bt which is one of the key values in 
this study. We are confirmed that the above four data sets are classified into two types 
according to each statistical property of the spread. Namely, the spread of USD/JPY 
or EUR/JPY exchange rates is time-dependent and fluctuates, whereas, the spread of 
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Fig. 1 Empirical distributions of 'return' (the difference between successive mid points) 
Amt = mt+i — rat for our four kinds of data sets. From the upper left to the 
lower right, P(Amt) for EUR/JPY exchange rates, USD/JPY exchange rates, Nas- 
daqlOO and price of gold are plotted. To compare the results with normal Gaus- 
sians, we calculate the empirical mean Am = (1/T) ^f-Tg Amt and the empirical 
variance = (1/T) JZ^Tg 1 (^Tit — Am) 2 for each data and plot the Gaussian 

m ) in the same panel. From the upper left to the lower right, these normal 
Gaussians are A/"(0.000012, 0.000018), A^(-0.000007, 0.000013), Af (-0.000726, 0.017722) and 
A/"(0.000080, 0.000922), respectively. 

NasdaqlOO or price of gold is a time- independent constant. In Fig. [2] for EUR/JPY 
and USD/JPY exchange rates, and in Fig. [3] for NasdaqlOO and price of gold, we plot 
the mid point nit, the return of the mid point Amt = "it+i — nit as a function of t 
and the distribution of the Bid- Ask spread P(S). From these figures, we clearly find 
that the Bid-Ask spread for the exchange rates apparently fluctuates, whereas the 
spread for NasdaqlOO or price of gold is a constant leading up to a single delta peak 
in the empirical distribution. From now on, the data in which the spread fluctuates 
is refereed to as data with stochastic Bid-Ask spread, whereas, the data in which the 
spread is constant is called as data with constant Bid-Ask spread. One of the main goals 
of this paper is to reveal the relationship between the statistical properties of Bid- Ask 
spread and the behaviour of auto-correlation and response functions for double-auction 
markets. 

3 Empirical data analysis 

In this section, we evaluate two macroscopic dynamical quantities, namely, auto-correlation 
and response functions by making use of empirical data analysis. These two relevant 
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Fig. 2 Statistical properties of the data set in which the spread fluctuates in time, EUR/JPY 
exchange rates (the upper panels) and USD/JPY exchange rates (the lower panels) are shown. 
From left to right, the mid point mt, the return of the mid point Amt = mt+i — mt as a 
function of t, and the distribution of the spread P(S) are plotted. 



xr 



PfSJ 




















Fig. 3 Statistical properties of the data set in which the spread is a time-independent constant, 
NasdaqlOO (the upper panels) and price of gold (the lower panels) are shown. From left to 
right, the mid point mt, the return of the mid point Amt = m t +i — mt as a function of t, and 
the distribution of the spread P(S) are plotted. 



quantities are explicitly denned by 

1 T_1 

C(l) = lim - ^t+l (1) 

t=0 

1 

R(l)= lim - Ve t (m t+; -mt). (2) 

f=0 

In order to evaluate these functions, we need the information about Selling-Buying 
signal et . However, as we mentioned in the previous section, the data gathered through 
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the MetaTrader4 |20] does not contain any information about it explicitly. To overcome 
this problem, we here assume that et is given in terms of 'return' of the mid point: 

et = sgn(m t+ i - m t ). (3) 

Namely, we assume when the mid point increases at the instant mj+i > mt, the 
number of traders who posted their own buying signal to the market also increases. As 
the result, the Selling-Buying signal e is more likely to take buying +1 at that instant 
t. 

Under the above assumption, for our four data sets, namely, EUR/JPY, USD/JPY 
exchange rates, NasdaqlOO and price of gold, we calculate the auto-correlation function 
C(l) and the response function R(l) via ([!} and ([2]), respectively. 



3.1 Auto-correlation function 

We first plot the auto-correlation function for the above four data sets in Fig. [4] From 




Fig. 4 Typical behaviour of the auto-correlation function C(l). From the upper left to the 
lower right, we plot EUR/JPY exchange rates (the estimated p = 0.24), USD/JPY exchange 
rates (the estimated p = 0.28), NasdaqlOO (the estimated p = 0.47), price of gold (the esti- 
mated p = 0.51). 

the upper left to the lower right, we plot EUR/JPY exchange rates (under the assump- 
tion on the asymptotic form: C{1) ~ p l , I > 1, the estimated p = 0.24), USD/JPY 
exchange rates (the estimated p — 0.28), NasdaqlOO (the estimated p = 0.47 ), price 
of gold (the estimated p = 0.51). From these panels, we find that the correlation in 
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the Selling-Buying signals decreases in the time difference I although the result for 
USD/JPY exchange rate possesses the negative correlation in I = 1 and converges to 
zero with a slight oscillation. 



3.2 Response function 

We next evaluate the response function for our four data sets. The results are shown 
in Fig. [5] In this figure, we plot the response function for the data with stochastic 




R(l) 



R(l) 




Fig. 5 Typical behaviour of the response function for our empirical data: from the upper 
left to the lower right, the results for NasdaqlOO, price of gold, EUR/JPY exchange rate and 
USD/JPY exchange rate are shown. 



Bid-Ask spread (the lower panels) and for the data with a constant Bid- Ask spread 
(the upper panels). From these panels, we find that some 'non-monotonic' behaviour 
in R(l) appears for the data set with stochastic Bid-Ask spread. 



3.3 Relationship between C(Z) and R(l) 

From the definitions, both auto-correlation and response functions are functions of the 
time-difference I. In the previous subsections, we investigated their behaviour inde- 
pendently. However, it might be assumed that these two quantities are related each 
other. Therefore, it is useful for us to make 'scatter plots' to reveal the dynamical 
relationship underlying these two quantities. In Fig. [6] we plot the relationship be- 
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Fig. 6 The relationship between R(l) and C(l). The upper two panels are results for the 
gold and NasdaqlOO with constant spreads, whereas the lower panel denotes the result for 
the EUR/JPY exchange rate having fluctuating spreads. We find that the linear relationship 
R(l) oc — C(l) holds for the data having a constant spread, whereas the linear relationship is 
apparently broken down in the EUR/JPY exchange rate which possesses a fluctuating Bid- Ask 
spread. 



tween R(l) and C(l) by scatter plots. The upper two panels are results for the gold and 
NasdaqlOO with constant spreads, whereas the lower panel denotes the result for the 
EUR/JPY exchange rate having fluctuating spread. We find that the linear relation- 
ship R(l) oc —C(l) is apparently broken down in the EUR/JPY exchange rate which 
possesses a fluctuating Bid-Ask spread. 

In the next section, we examine a phenomenological model to explain the non-linear 
relationship R(l) oc —f{C(l)) (f{x) denotes a non-linear function) theoretically. 



4 A phenomenological approach 

In order to explain the behaviour of the auto-correlation and response functions, we 
examine a phenomenological approach based on the so-called Madhav an- Richards on- 
Roomans model (MRR for short) |15II16| to simulate the stochastic process of price- 
change in empirical data sets. In the MRR model, the price pt updates according to 
the following rule. 

Pt+l = Pt + 0(et - pet-l) + & (4) 

where £t denotes a noise in the market satisfying = and (£t£ t '} = 8 t t '^ 2 - 6 is a 
constant value to control the slope of instantaneous price change. The label et means a 
Selling-Buying signal to represent et = +1 for a buying signal and et = — 1 vice versa. 
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Behaviour of the above update rule is dependent on the statistical properties of 
Selling-Buying signals et- p is a correlation factor and in the MRR theory, we assume 
that et follows a simple Markovian process, namely, 

J2 e t P(et\e t -i) = pet-v (5) 

e t =±l 



The price value of pt+i provided that the Selling-Buying signal in the previous 
time step is et = +1 should be the Ask at and the price pt+i provided that the signal 
is et = — 1 should be the Bid bt- Hence, we naturally define the time-dependence of 
Ask and Bid as follows. 

at = p t + 9(1 - pe t -i) + <j> (6) 
b t = pt + 9(-l - pet-i) - 4> (7) 

where (j> denotes a kind of transaction cost and the value itself is set to a constant in 
the MRR model. From these rules, we easily find the Bid-Ask spread at time t as 

S t = a t - 6* = 2((9 + </>). (8) 

Namely, in the MRR model, the spread is a time-independent constant during the 
dynamics. 

On the other hand, the mid point of the Bid and Ask is given by 

mt = ^(at + h) =pt~ 9pe t -i- (9) 

Therefore, for the parameter choice 6 — or p — 0, the mid point mt is identical to 
the price pt- For the above update rules of price, Bid, Ask, spread and mid point, we 
investigate the macroscopic properties of double-auction markets through the auto- 
correlation function and the response function. 



4.1 Auto-correlation function 

From the definition of Markovian process ((SJ), the auto-correlation function is given by 

C(l) = {etet+i) = Yl "■ J2 P^t,--- ,e t+ i)ete t+ i^p l . (10) 
£ t =±l e t+! =±l 

We should keep in mind that the auto-correlation function is originally defined by 
JTJ). However, in the limit of T — > oo, one can replace the time-average by the average 
over the joint probability of the stochastic variables et, • • ■ , e t +i as (1/T) Y1I'=q(' ") = 
J2 et =±i ■ ■ ■ J2e t+l =±i(- ■ ■ ) P ( e «> ' ' ' i e t+l) according to the law of large number. 

The correlation factor p should be \p\ < 1. For a positive p, the correlation function 
decays exponentially as C(l) = e ~ n ° s(1/p K 
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4.2 Response function 

We next consider the response function of the market, that is defined by 

R(l) = (et ■ (m t+ ; - m t )) (11) 

where the bracket (• ■ ■ } has the same meaning as that in (|10l) has. 

From the above response function, one obtains some information about the response 
of the market at time t + I to the Selling-Buying signal at arbitrary time t. Namely, 
the response function measures to what extent the mid point increases (decreases) 
on average for interval I when a buying (selling) signal is posted to the market I steps 
before we observe the mid point. After simple algebra, we easily obtain the relationship 
between the response function R(l) and the correlation function C(l) for the MRR 
model as follows. 

R(l) = 8(1 - C(l)) (12) 

As we saw before, for a positive correlation factor p > 0, the C(l) monotonically de- 
creases as C(l) = e - a °s(Vp) . Hence, the response function also behaves monotonically 
and converges to 8 as R(l) = 6>(1 - e~ Hog(1/p) ) ->• 8 (I -> oo). It should be noticed that 
the linear relationship between C(l) and R(l) holds from (|12ll . 

From equation (|12p . we also find -R(l) = 8(1 — p) and this fact tells us that 
R(oo)/R(l) = (1 - p) -1 holds. The volatility defined by 

a 2 (l)^j((m t+l -m t ) 2 ) (13) 

also reads 

a\l) = E 2 + e\l - pf jl + | (14) 

and a 2 (l) = E 2 + 8 2 (1 - p) 2 , a 2 (oo) = E 2 + 8 2 (1 - p 2 ). We may use the above rigorous 
equations to check the validity of our computer simulations. 

In Fig. [7] we show the typical behaviour of response function for several choices 
of p. From these panels, we find that for a positive correlation factor, the response 
function monotonically converges to the value 8. In this figure, we also show the results 
obtained by simulating the update equation for the price ([4| and the mid point ((9]), 
and calculating the response function numerically by making use of ©. We find that 
the both theoretical prediction (solid lines) and the simulation (boxes) are in good 
agreement. Moreover, for the choice of positive correlation factors, say, p = 0.1 and 
0.5, the response function increases monotonically and converges to 8(— 1.5) as the 
MRR theory predicted. 

In Fig. [8] we show the relationship between R(l) and C(l) for the MRR model. We 
clearly find that the linear relationship R(l) = #(1 — C(l)) oc — C(Z) actually holds. We 
should notice that this result is completely different from the relationship shown in the 
lower panel of Fig. [()] which is result for the data having fluctuating Bid- Ask spread. 
This failure of the MRR model to simulate the 'A-shape' of the C(l)-R(l) scatter plot 
for the data with stochastic Bid-Ask spread is obviously due to the assumption of 
constant Bid- Ask spread, namely, St = at — h = 2(8 + <f>) on the MRR theory. We also 
conclude that the breaking of the linear relationship R(l) oc — C(Z) is macroscopically 
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Fig. 7 Typical behaviour of the response function for the MRR model. From the upper left 
to the lower right, we plot the response function for p = —0.9, —0.1, 0.1, 0.5. The solid lines are 
theoretical predictions by R(l) = 9(1 — p )■ The boxes are obtained by numerical simulations 
for a finite T = 10 5 . 

due to the fact that the response function evaluated on the basis of the MRR theory 
behaves 'monotonically' and converges to a finite value 8. 

This fact is one of the limitations of the MRR theory to explain the empirical 
data for double-auction markets. To make the efficient model for the double-auction 
market with stochastic Bid- Ask spread, we use a game theoretical approach based on 
the so-called minority game. 



5 A minority game modeling of double-auction markets 

In order to make a model to simulate the 'non-monotonic' behaviour of the response 
function for the financial data with stochastic Bid-Ask spread, we start our argument 
from standard minority game 22 , 23 , 24 , 25 with a finite market history length. 



5.1 General set-up 

In our computer simulations for the minority game, at each round t (time step) of 
the game, each trader i (i = 1, • • • , N: N should be an odd number to determine the 
'minority group') decides his (or her) decision: B{(t) = +1 (buy) or Bi(t) — —1 (sell) 
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Fig. 8 The relationship between R(l) and C(l) for the MRR model. We find that the linear 
relationship R(l) = 9(1 - C{1)) <x -C(l) holds. 



to choose the minority group. Then, we evaluate the total decision of the traders: 



N 



A(t) 



(15) 



i=i 



for each round t. It should be noted that the factor iV -1 / 2 is needed to make the A(t) 
of order 1 object (independent of the size N). From the definition (|15[) . the market 
is seller's market for A(t) > 0, whereas the market behaves as buyer's market for 
A(t) < 0. The A(t) follows complicated stochastic process and one might consider that 
the price p is updated in terms of the A(t) as follows. 



p{t + 1) = P (t) + 0{A(t) + V sgn[A(t - 1)]} 



(16) 



where j3 and V are positive constants. The above update rule means that the price 
increases if the 'buying group' is majority and decreases vice versa. This setting of 
the game seems to be naturally accepted. A bias term ip appearing in (|16p plays an 
important role to simulate the auto-correlation function as we will see later on. 

To decide 'buy' or 'sell', each trader uses the following information vector defined 

by 



X(A, Z : t) 



( sgn[(l-C)i4(t-l)+CZ(t,l)] \ 
sgn[(l - t)A(t - 2) + {Z(t, 2)] 



(17) 



where sgn(x) denotes a sign function and Z(t, £), £ = 1, • ■ 
by 

(z(t,0z(t',O) = s tt ,s u ,. 



\ Sg a[(l-C)A(t-M)+CZ(t,M)]J 

, M is a white noise defined 



(18) 



Therefore, each trader uses the information of market through the up-down configura- 
tion of the return A (with some additive noise Z) back to the previous M-steps. If the 
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£ is close to 1, the 'real' market history through the A is hided by the 'fake' market 
history through the noise Z. 

For a given information vector X(A, Z : t) — A chosen from all possible 2 M can- 
didates, each trader i decides her (or his) action at round I by the following strategy 
vector: 



i i il 



,r*A), A = l, 



(19) 



where we defined A as the index (entry) of the selected information vector A and s 
stands for the number of the possible strategies for each trader. Each component of 
the above strategy vector r\ takes +1 (buy) or —1 (sell). Therefore, each trader has 



her (his) own look-up table which is defined by a matrix with size s x 2 as 



R 



r 2 



rf 
rf 



r-i 



\ 



(20) 



Each component of the above look-up table = — 1} is fixed ('quenched') before 
playing the game. However, in the next section, we consider the case in which the com- 
ponents of look-up tables are rewritten during the game. In this paper, we concentrate 
ourselves to the simplest case of two strategies s = 2 and £ = ('real' market history). 
Then, the trader i changes her/his own pay-off value Vi c by the following update rule: 



Vi c (t + 1) = V ic (t) 



Bi(t)A(t) 



Bi{t) = 8(A,X(A,Z: t))r 



ici{t) 
A 



(21) 
(22) 



A=l 



where S(x,y) denotes the Kronecker's delta and Cj(t) means the optimal strategy in the 
sense that Cj(t) is given by 



Ci(t) = argmax c [P ic (t)]. 



(23) 



The meaning of the update rule ()21|) is given as follows. If A(t) > and the majority 
group consists of traders who post their decisions +1 to the market, the trader i 
attempts to post her/his decision as an opposite sign of A(t), namely, Bi(t) = — 1. 
Thus, the trader i acts so as to satisfy the condition Bi(t)A(t) < which leads to 
increase of her/his pay-off value Vi c (t + 1). 

By taking into account the fact that we are dealing with the case of s = 2 (c = 1, 2), 
we rewrite the equation (|2ip as 



g t (t + 1) = qi (t) - -L V S(A, X(A, Z : t))rf A A{t) 

VN A=l 



(24) 



by means of 



(25) 
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Substituting (1221) into the definition of the total bit A(t), we have 
N N 2 M 

i=l i=l yl=l 

We should notice that the above equation can be written by using the following relation 



r~* " (t) =t"A + sga[qi(t)]vA, u>A = \{ r A + (27) 



Then, we obtain the following coupled non-linear equations with respect to the total 
decision A(t), the difference of pay-off values for two strategies <ft(t) and the update 
equation of the price p(t): 



iv z 

A ^ = ^^^S(A,X(A,Z:t))[w A + S ga[ qi (t)]r, A \ (28) 

^ i=l A=l 

©(« + 1 ) = 3iW-4E ^ A ( A ' Z : ( 29 ) 
VJV A=l 

P (t + 1) = p(t) + + V sgn[A(t - 1)]}. (30) 

The above rules (|29I)(I28I) and (|30|) are our basic equations to discuss the response 
of double-auction markets having stochastic Bid-Ask spread to instantaneous Selling- 
Buying signals. 



5.2 Making of the Bid-Ask spread in the minority game 

To make the Bid-Ask spread in our minority game, we assume that the buying price 
an and the selling price 6j t which are posted to the market by each trader i at round 
(time) t are updated according to the following rules. 

"it = Pif) + la gu + S (31) 
kt = p(t) + lb 9it - 5 (32) 

where 7a,7fc and S are constants to be set so as to satisfy an — bn > 0. gn is an 
uncorrelated Gaussian variable with mean {git) = and covariance {gags *') = 5 t t i& i ./ 
(additive white Gaussian noise: AWGN). In our simulations, we set ~y a = lb — 0.01, 5 — 
0.049. Then, the Bid-Ask spread at round t is given by 

S t = min{a lt |a it £ N+} - max{b lt j6 it G 7V_} (33) 

where the groups taking 'buying' and 'selling' decisions are refereed to as N+ and N-, 
respectively (N = N+ + N-). 

In Fig. [9] we plot the resulting distribution of the spread S — St- From Fig. [9l we 
find that the stochastic Bid- Ask spread generated from the above modeling based on the 
minority game actually fluctuates and possesses a non-trivial shape of the distribution. 
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0.009 




0.055 



Fig. 9 Distribution of the spread S = St generated by the minority game. 
5.3 Results 

For the above set-up of the minority game, we evaluate two relevant statistics, namely, 
correlation function C(t) (by ([1])) and the response function R(t) (by p])) to compare 
the results with the empirical evidence for the data with stochastic Bid- Ask spread. 

5.3.1 Auto-correlation function 

We first examine the effect of the bias term ip on the correlation function. The results 
for if) — are shown in Fig. [TU] (left). From this panel, we find that C(l) — for I > 2 




2 3 4 5.6 7 8 9 1 2 3 4 5.E 7 3 9 



Fig. 10 Auto-correlation function C(l) in the minority game for ip = (left). We set TV = 
1025, M = 9, ft = 0.01 and iterated the game 100010 rounds. The C(l) is zero for I > 2. The 
right panel is auto-correlation function of the MRR model with p = 0.1. The error-bars were 
calculated by 10-independent trials. 



and the result is apparently different from the result for the empirical data. Here we 
set N = 1025, M = 9, /3 = 0.01 and iterated the game t = 100010 rounds. The right 
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panel of Fig. 1101 shows the correlation function of the MRR model for p — 0.1. The 
Selling-Buying signal 



e t = sga[A(t)] 



(34) 



is actually correlated automatically through the market history with length M = 9, 
however, the correlation strength is very weak. Therefore, we need some other explicit 
correlation through the bias term ip which enhances the correlation by two-round back 
A(t — 1) from the present 4+1. Hence, we here choose the non-zero bias term ip to 
reproduce the auto-correlation function as observed in the empirical data. 

We checked that the results are robust against the slight differences in the param- 
eters appearing in the game such as /3,M,S etc. However, for only parameter ip, we 
should be careful to choose the value. This is because from the definition of update 
rule of the price (|30[) . the effect of the A(t) on the price change is relatively depressed 
by the large value of the ip. Therefore, we should choose ip so as to make the value 
smaller than the standard deviation of the A (t), namely, square root of the volatility 



ip < 



A(t)Y 



a A . 



(35) 



In our simulation, the square root of the volatility is estimated as a a = 0.44295. In Fig. 
1111 we plot the correlation function for the case of ip — 0.05 (left panel) and ip — 0.1 
(right). We should notice that these two choices of the bias term ip satisfy the condition 
J35J. 




Fig. 11 Auto-correlation function C(l) in the minority game for ip = 0.05 (left) and ip = 0.1 
(right). We set N = 1025, M = 9,/3 = 0.01 and iterated the game for 100010 rounds. The 
error-bars were calculated by 10-independent trials. 



From these panels, we find that the correlation function decreases as we observed 
in the same function for the empirical data sets. 

5.3.2 Response function 

We next plot the response function in Fig. [T^] for ip — 0. (upper panel), ip — 0.05 
(lower left) and ip — 0.1 (lower right). From these panels, we find that the behaviour 
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R(l) 




Fig. 12 Response function R(l) in the minority game for tp = (upper), ip = 0.05 (lower 
left) and if> = 0.1 (lower right). We set N = 1025, M = 9, /S = 0.01 and iterated the game for 
100010 rounds. The error-bars were calculated by 10-independent trials. 



of the response function is not monotonically increasing function leading up to the 
convergence to some constant value but 'non-monotonic' as the response function of 
data sets having stochastic Bid- Ask spread (EUR/JPY, USD/JPY exchange rates) 
shows (see the lower panels of Fig. O . 

5.3.3 Relationship between the auto-correlation and response functions 

In Fig. 1131 we plot the relationship between the auto-correlation and response functions. 
In the upper panel, we show the result for zero bias term tjj — 0. we find that the curve 
is deviated from the linear relation R(l) oc — C(l). However, the shape is not 'A' as 
observed in the empirical evidence but T-shape'. In the lower two panels are results 
for the non-zero bias term ijj ^ 0. We clearly find that the 'A-shapes' appears and the 
results are qualitatively similar to those of the empirical data. 



6 Adaptive look-up tables 

In the previous section, we fixed each decision component (buying: +1, selling: —1) in 
their look-up tables before playing the game (in this sense, the decision components are 
'quenched variables' in the literature of disordered spin systems such as spin glasses). 
However, in this paper, a certain amount of traders update their decision components 
according to the macroscopic market history (they 'learn' from the behaviour of mar- 
kets) so as to be categorized into two-groups with a finite probability (in this sense, 
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Fig. 13 The relationship between CO.) and RO) for minority game with ip = (upper), 
if, = 0.05 (lower left) and if) = 0.1 (lower right). We set N = 1025, M = 9,/3 = 0.01 and 
iterated the game for 100010 rounds and we calculate C(i),i?(Z) up to I = 100. 



the components are now regarded as 'annealed variables'). Namely, at each round of 
the game, if the number of sellers is smaller/greater than that of buyers, a fraction 
of traders, what we call optimistic group / pessimistic group, is more likely to rewrite 
their own decision components from — 1/+1 to +1/— 1. In realistic trading, we might 
change our mind and rewrite the components of the look-up table according to the 
market history. Therefore, in this section we consider the case in which some amount 
of traders can rewrite their own table adaptively. 



6.1 Adaptation using the latest market information 

We first consider the case in which each trader updates her/his own look-up table 
according to the latest information of the market. Some of the traders make their 
decisions as 'buy' when the market is seller's market, namely, the signal of the latest 
market is 'buying' (what we call optimistic group). On the other hand, they decide 'sell' 
vice versa (they are referred to as pessimistic group). Namely, each trader rewrites the 
table according to the following algorithm. 

Adaptation algorithm using the latest market information 

(i) Fix ('quench') each component of the look-up table at the beginning of the game 
t = -M. 
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(ii) At each game round t for t > —M, each trader rewrites her/his component r l \,r % A , 
where A denotes the entry of market history for the information vector \(A, Z : t), 
with probability fa as 

r 1 } = Bgn[(A(t - 1)], = sgn[A(* - 1)]. 

(iii) Each trader recovers her/his original (at the beginning of the game) look-up table 
with a probability fa at the next game round t + 1. 

(iv) Repeat (ii) and (iii) until the game is over. 

Namely, a fraction ~ Nfa of the traders is categorized into the 'optimistic group' if 
A(t — 1) > (seller's market) and into the 'pessimistic group' if A(t — 1) < (buyer's 
market). 

6.1.1 Results 

In Fig. 1141 we show dynamical quantities C(l) (upper left, middle left) and R(l) (upper 
right, middle right) evaluated for the minority game with adaptive look-up tables. We 
set ip = 0, fa = 0.01 and fa = 1 (upper panels), fa — 0.9 (middle panels). The lower 
two panels show the relationship between the auto-correlation and response functions 
for (fi,fa) = (0.01, 1) (left) and (fa, fa) = (0.01, 0.9) (right). From this figure, we find 
that the positive correlation C(l) > for I — 1 appears even if we set ip = and non- 
monotonic behaviour of the response function is reproduced. As the result, 'A-shape' 
in the C(l)-R(l) scatter plots is generated. From these results, we conclude that the 
adaptive modification of the look-up table by using the latest market information for 
each trader works well to explain the empirical evidence. 



6.2 Adaptation by using the market history 

In the previous subsection, we succeeded in generating a positive finite auto-correlation 
by making use of the adaptive look-up table even if we set ip — 0. However, in this look- 
up table, each trader changes her/his decision from the latest information about the 
market. As the result, the auto-correlation function decays to zero for / > 2. To modify 
the weak correlation, we construct the adaptive look-up table by using the information 
about the market history with length M. As we mentioned, the information vector 
X(A, Z : t) contains the useful information on the market. Hence, we shall assume that 
each trader rewrites the component of her/his table as 



with probability 



r l l =sga[n A (A,Z:t)], r%=sgn[f2 A (A,Z:t)] (36) 



a\n A (A,Z:t)\ 
« = ( 37 ) 



M 



where we defined 

Q A (A, Z :t) = w- \(A, Z :t) = Y^ 2 M ~ T sgn{A(t - r)] (38) 

T = l 

w = (2 M-1 ,2 M ~ 2 ,--- ,2°) (39) 
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Fig. 14 Dynamical quantities C(i) (upper left, middle left) and R(l) (upper right, middle 
right) evaluated for the minority game with adaptive look-up tables. We set ip = 0, /i = 0.01. 
fs = 1 (upper panels), = 0.9 (middle panels). The lower two panels show the relationship 
between the auto-correlation and response functions for (/i, $2) = (0.01, 1) (left) and (/1, f2) = 
(0.01,0.9) (right). 



at each game round t. We should keep in mind that each trader recovers her/his original 
look-up table with a probability fa at the next game round. The i?yi(^4, Z : i) denotes 
cumulative weighted market status, namely, we assume that the importance of the 
market information decays as 2 _T in the history length r. 

For instance, for the information vector (let us define the entry by A) having +1 for 
all components: X(A, Z : t) = (sgnL4(t - 1)], • • • , sgn[A(t - M)]) = (+1, ■ ■ ■ , +1), that 
is, the market remains as seller's market up to M-times back, we obtain Qj\{A,Z : 
t) = 2 M — 1 > and iV/3 = iVa(l — 2~ M )-traders rewrite their components as 
r A' r A = s g n [^/l(A Z, t)] = 1. As another example, when seller's market and buyer's 
market appears periodically as \(A, Z : t) = (—1, +1, —1, +1, ■ ■ ■ ), we have f2/i(A, Z : 
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t) = -{2 M - (-l) M /2}/3 < and Nfy = Na{l - (-l) M /2 M+1 }/3-traders rewrite 
their component as r A ,r A = sgn[f2 A (A, Z : t)] = — 1. 

Let us summarize the above procedure as the following algorithm. 

Adaptation algorithm using the market history 

(i) Fix ('quench') each component of the look-up table at the beginning of the game 
t = —M. 

(ii) At each game round t for t > —M, each trader rewrites her/his component r A ,r A , 
where A denotes the entry of market history for the information vector X(A, Z :t), 
with probability as 

rl = sgn[0 A (A,Z: t)], rf = sgn[!2 A (A, Z, t)] 

with probability 

_ a\n A (A,Z:t)\ 
J 3 2 M 

with Q A (A, Z:t) = J2r=l 2 M ~ T sgn[A(t - r)}. 

(iii) Each trader recovers her/his original (at the beginning of the game) look-up table 
with a probability fa at the next game round t + 1. 

(iv) Repeat (ii) and (iii) until the game is over. 




We set a = 0.01 and select the same values as those of the previous section for the 
other parameters. 

6.2.1 Results 

In Fig. 1151 we first plot the generated distribution P(f2) of f2 A (A,Z : t) = Q (left) 
and the typical time-evolution of the probability fa for the first 200 steps (right). It 
should be noted that Q in the left panel is ranged from i? m ; n = X^r=l 2 M ~ T (— 1) = 
1 - 2 9 = -511 to fl mas = Y!^=\ 2 M ~ t (+1) = 2 9 - 1 = 511 because we set the history 
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Fig. 16 Dynamical quantities C'(l) (upper left, middle left) and R(l) (upper right, middle 
right) evaluated for the minority game with adaptive look-up tables for M market history 
length. We set i)> = 0, fa = 1 (upper panels), f 2 = 0.9 (middle panels). The lower two panels 
show the relationship between the auto-correlation and response functions for fa = 1 (left) 
and fa = 0.9 (right). 



length M — 9 in the definition of the 0^(A, Z : t). We find from the right panel that 
the typical time-evolution of the probability f$ obeys complicated dynamics. 

We next show the results for the macroscopic quantities in Fig. 1161 From this 
figure, we confirm that the 'A-shape' in R(l)-C(l) scatter plots are much similar to the 
empirical evidence than the results in the previous subsection. 



7 Discussion 



Our general set-up presented in this paper is applicable to the analysis for the other 
quantities or for the other stochastic models. Here we shall mention them briefly. 
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7.1 Waiting time statistics 

As we mentioned in section [T] we can generate the duration between the price changes 
within the framework of our minority game. 

Let us introduce the maximum value of the spread S which is determined by market 
makers. Then, we might define a set of time points at which the price is updated as 

{t} = {t k \S tk = mm{a ltk \a ltk g N+} - max{b itk \b itk e N-} < S} (40) 

For these time points t k , the duration between successive price changes is given by 

M = {r k \r k = t k+1 - t k } (41) 

We can evaluate the distribution P(t) and compare the results with well-known dis- 
tributions, for instance, the Mittag-Leffler type 26,27,28 . 



7.2 Mean-field models 

Recently, Vikram and Sinha 30 proposed a mean-field model to describe the collective 
behaviour of financial markets. In their model, each trader i decides his (her) decision: 
Si(t) = +1 (buy), —1 (sell) and Si(t) = (no action) at time t according to the 
following probability. 

P[\Si(t)\ = 1] = 1 - P[Si(t) = 0] = exp (-n 

where p(t) stands for the price at time t and the bracket (• ■ • ) T means the moving 
average over the past r-time steps. The parameter fj, is a parameter which controls 
the sensitivity of an agent to the magnitude of deviation of the price from its moving 
average (p(t)). For fi = 0, the system reduces to a binary decision model where every 
agent trades at all time instants. Namely, P[Si(t) = 0] = for fi = means that the 
trader decides his (her) action, buy or sell, certainly. The traders who decide to trade 
at t make his (her) action randomly, that is, P[Si(t) = 1] = P[Si = —1] = 1/2. The 
price at time t + 1 is decided by the following recursion relation. 

Pit+l)= {\ ± ^l) p{t) (43) 

where At = (1/s/N) y~)j?-\ Si(t). In Fig |17l we plot the typical time-evolution of the 
price p(t) and the moving average (p(t)) T evaluated by (|42"T) and (|43)l . 

To construct the double-auction market for the above price change, we shall use 
the same definitions of buying bn and selling an signals for each trader i at time t as 
pip (|32[l in our minority game modeling. Then, we calculate the response function and 
the auto-correlation function by using the above set-up and plot them in Fig. 1181 We 
set the number of agents JV = 20000, the number of iterations for the price change 
T = 100000. We also choose 8 = 0.02, 7 a = j b = 0.001. The initial condition on the 
price p(t) is chosen from the Eur/JPY exchange rate in the empirical data. From this 
figure, we find that the behaviour is different from the empirical evidence. Especially, 
the auto-correlation function fluctuates with very large amplitudes. This result might 



log 



P(t) 
(P(t)h 



(42) 
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Fig. 17 Typical behaviour of time evolution of the price p(t) (the solid line) and the moving 
average {p(t)) T (the broken line) evaluated by ll42l and ll43l l. We set fi = 100, r = 10000. 



Fig. 18 Typical behaviour of the response function (upper panels) and the auto-correlation 
function (lower panels) for the Vikam-Sinha model. From the left to the right, we set (fi, t) = 
(100, 10000), (10, 10000), (100, 1000). The error-bars were calculated by 10-independent trials. 



come from the fact that in the Vikam-Sinha model, the buying-selling signal is chosen 
randomly. Therefore, we should modify the Vikram-Sinha model in order to explain 
the non-monotonic behaviour of the response function in double-auction markets with 
stochastic Bid-Ask spread. 



8 Summary 

In this paper, statistical properties of double-auction markets with Bid- Ask spread were 
investigated through the response function. We first attempted to utilize the so-called 
Madhavan-Richardson-Roomans model (MRR for short) to simulate the stochastic pro- 
cess of the price-change in empirical data sets (say, EUR/JPY or USD/JPY exchange 
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rates) in which the Bid- Ask spread fluctuates in time. We found that the MRR theory 
apparently does not simulate so much as the qualitative behaviour ('non-monotonic' 
behaviour) of the response function calculated from the data. Especially, we were con- 
firmed that the stochastic nature of the Bid-Ask spread causes apparent deviations 
from a linear relationship between the R(l) and the auto-correlation function C(l), 
namely, R(l) oc —C(l). To make the microscopic model of double-auction markets 
having stochastic Bid- Ask spread, we utilized the minority game with a finite market 
history length and found numerically that appropriate extension of the game shows 
quite similar behaviour of the response function to the empirical evidence. We also 
revealed that the minority game modeling with the adaptive ('annealed') look-up table 
reproduces the non-linear relationship R(l) oc —f(C(l)) (f(x) stands for a non-linear 
function leading to 'A-shapes') more effectively than fixed ('quenched') look-up table 
does. 

Of course, there are still gaps between the theoretical prediction and the empirical 
evidence. We should modify our modeling of the double-auction market and figure out 
the micro-macro relationship in the market much more quantitatively. 
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